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Abstract. The interplay between stellar populations and gas in local starburst galaxies is 
analyzed using images from the Hubble Space Telescope to map the ages of the young stellar 
components and to isolate the contribution of shocks on spatial scales ranging from a few tens 
of pc to ~1 kpc. The shocked gas represents a small fraction of the total ionized gas in these 
objects, yet it can have profound effects on the long-term evolution of the starburst, which may 
include the triggering of new star formation. 



1. Introduction 

The evolution of galaxies and of the intergalactic/intracluster medium (IGM/ICM) are 
closely connected, via the continuous interchange between the ISM and the IGM/ICM. 
The heating and metal enrichment of the IGM/ICM, and the enrichment of galaxies' 
ISM with pristine gas via infall from the IGM are manifestations of this close con- 
nection. Evidence for gas outflows from galaxies has been found at redshift as high as 
~3 l|Pettini et al. 20 001. and may have contributed to the early pollution of the IGM 
IjEllison et al. 20000 . 

The driving mechanism for the energy and metal pollution of the IGM/ICM is the 
feedback from star formation. In the basic scenario, the combined effects of ejecta and 
energy release from supernovae and massive stars in regions of intense star formation form 
super bubbles that sweep up the surrounding medium and produce cavities containing hot, 
shocked gas (Weaver et al. 1977). If the bubble has enough energy to expand to a size of 
a few times the scale height of the ambient medium, the swept-up shell may break out 
and the wind internal to the bubble may expand into the IGM (Chevalier & Clegg 1985). 

Many theoretical and observational studies have investigated the mechanical feedback 
from star formation and its efficiency in transfering energy and material to regions far 
removed from the site of star formation (e.g., De Young & Heckman 1994; Suchkov et al. 
1994; Mac Low & Ferrara 1999; Strickland & Stevens 2000; Heckman, Armus & Miley 
1990; Dahlem 1997; Ferrara & Tolstoy 2000; Cecil, Bland-Hawthorn & Veilleux 2002), 
as well as in triggering renewed star formation (e.g., Wada & Norman 2001). Yet many 
questions remain to be fully answered. Among these, progressing from small to large 
scales: the impact of the spatial distribution, history, and duration of the star-forming 
event (Clarke & Oey 2002); the coupling of the feedback from star formation with the 
surrounding ISM; the dependence of the feedback efficiency on the conditions of the ISM 
(e.g. porosity; Clarke & Oey 2002); and the role of the global galaxy parameters (mass, 
metallicity, environment; Dahlem 1997; Strickland et al. 2004). A review of some of these 
issues is given by D. Strickland (2004, these proceedings). 

One of the complicating factors in tackling the various facets of feedback is the vast 
range of physical scales that needs to be probed. While the relevant scales for charting 
star formation are those of the stars and star clusters (a few pc) , a study of the structure 
and conditions of the ISM needs to sample scales of tens to hundreds of pc. Moving on to 
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NGC3077 


3.85 


1.4 


-17.5 


8.9 


0.076 


0.043 


0.51-1.40 


NGC4214 


2.94 


2.0 


-17.2 


8.2 


0.089 


0.041 


0.29-1.24 


NGC5236 


4.5 


0.75 


-20.3 


9.2 


0.308 


0.030 


1.27-2.90 


NGC5253 


4.0 


1.5 


-17.5 


8.2 


0.270 


0.032 


0.38-1.19 



Table 1. Properties of the observed galaxies. Columns are as follows; D: distance; d: 
physical size of the region observed; Mb- absolute B magnitude; (O/H): oxygen abun- 
dance, 12+log(0/H); SFR(Ha): star formation rate derived from the extinction corrected 
Ha luminosity; L(Ha) np h/L(Ha): fraction of the Ha luminosity associated with shocks; 
L(Ha)^ p e h I L(Ha) np h'- ratio of the predicted-to-observed Ha luminosity associated to shocks. 
NGC3077 is in the M81 group, in close interaction with M81 itself; NGC4214 is an isolated 
galaxy, while NGC5236 and NGC5253 form a loose pair. 



larger scales, superwinds are best investigated on kpc scales, while probing the influence 
of global parameters involves scales of tens of kpc (galactic sizes) to Mpc (interactions) . 

Even for the closest galaxies, the small-to-intermediate scales relevant for the interac- 
tion feedback-ISM are accessible only via high angular resolution observations (e.g., 
with the Hubble Space Telescope), as 10 pc=0.4" for a galaxy at 5 Mpc distance. 
The study reviewed here involves four local starburst galaxies observed with the Wide 
Field Planetary Camera 2 on the HST. Because of their intense-to-extreme star for- 
mation rates, starburst galaxies are the sites where stellar feedback can have its most 
dramatic influence on the structure and evolution of the ISM (and surrounding IGM; 
|Heckman, Armus fc Miley 19901 ), and are, therefore, optimal laboratories to investigate 
mechanical feedback. The four galaxies in the present study cover a range in luminosity, 
metallicity, star formation rate, and environment, but they are all closer than 5 Mpc (Ta- 
ble^l. This investigation attempts at addressing the issue of the coupling feedback-ISM, 
and its relation to the duration of the star forming event. The problem is tackled from 
two fronts: stellar populations and ISM conditions. The ages of the young (<300 Myr) 
stellar clusters are used to set a minimum value to the duration of the current starburst 
event, while measurements of the shocked gas provide a constraint on the fraction of star- 
burst mechanical output that is recovered in the ISM within a small distance (<1 kpc) 
of the site of star formation. 

2. Charting the Young Stellar Populations in Starbursts 

The star cluster populations and the spatial distributions of their ages within the 
starburst sites were investigated for three of the galaxies listed in Table ^ NGC3077, 

NGC5236, and NGC5253 l|Harris et al. 20011 IHams et al. 2004|). For NGC5253, the dif- 

fuse starburst population had previously been age-dated IjCalzetti et al. 1997llTremonti et al. 2001|l . 
which, combined with the more recent study of its cluster population, yields a complete 
picture of the recent star formation history in this galaxy. 

Images of the galaxies in the HST/WFPC2 medium and broad filters UV, V, I, and 
in the Ha and H/3 line emission provide enough angular resolution to identify clusters, 
and enough color information to age-date them. In particular, two age-sensitive diag- 
nostics, the color-color diagram UV— V versus V— I and the equivalent width of Ha, 
were combined to constrain the most likely age for each stellar cluster and, assuming a 
Salpeter IMF, to assign a mass f Figure HI [Harris et al. 20011 IHarris et al. 2004jl . Ha/H/3 
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ratio maps were used to correct photometry and fluxes for the effects of dust redden- 
ing (Cal zetti, Kinney fc Storchi-Bergmann 1994| ). The cluster ages were then re- mapped 
onto the galaxy images to derive the history of the cluster formation in each starburst. 

The two dwarf galaxies, NGC3077 and NGC5253, show time-extended star formation, 
spanning the last 100-300 Myr, and a very 'chaotic' spatial distribution of the cluster 
ages. This similarity between the two galaxies is counterbalanced by striking differences. 
In NGC3077, there are star clusters as old as 300 Myr, while in NGC5253 the oldest 
clusters unambiguously detected in the HST images do not seem to be significantly older 
than ^20 Myr IjHarris et al. 2 004). Indeed, the time-extended star formation history 
for NGC5253 is inferred not from its clusters, but from its diffuse stellar population, 
whose colors are consistent with constant star formation over the past ^100-200 Myr 
IjCalzetti et al. 1997}. The dearth of 'old' star clusters in NGC5253 could be caused by 
rapid cluster dissolution timescales in the center of the galaxy; by rescaling the Milky Way 
model of |Kim, Morris &T Lcc 1999] lTremonti et al. 200ll and lHarris et al. 2004 derive dis- 
solution timescales in the range 16-50 Myr for NGC5253. These short timescales are 
driven by the high central velocity dispersion measured in the galaxy ( |Caldwell <k Phillips 1989 
A smaller velocity dispersion may be present in NGC3077, as inferred from molecular 
cloud velocities IjWalter et al. 2002)l . implying longer survival times for its clusters. 

Contrary to the dwarfs, the distribution of the stellar clusters in the center of the 
giant spiral NGC5236 shows a high level of organization, possibly a reflection of the 
different dynamical environment (presence of a bar) in the massive galaxy. The UV-bright 
clusters are distributed along a half-ringlet surrounding the optically-bright nucleus. The 
distribution of the ages also shows a high level of correlation in the star formation, which 
appears propagating in the S-N direction and inside-out in the ringlet (Harris et al. 
2001, |Puxley, Doyon fc Ward 1997| ). The re is a clear peak in the number of clusters in 
the narrow age range 5-7 Myr, with a handful of older clusters (up to ^50 Myr). Although 
the diffuse population is generally consistent with constant star formation over the past 
~1 Gyr, the sharp cut-off of the cluster ages around 7 Myr could still be explained with 
either a recent starburst or very rapid cluster dissolution timescales IjHarris et al. 200l)l . 

Independently of whether the current starburst in the grand design spiral is confirmed 
to be younger than ~10 Myr, a clear mechanism for the triggering and feeding of the 
starburst, i.e. gas infall along the main bar, is present in NGC5236 (Gal lais et al. 199l)l . 
No such obvious mechanism is present in the two dwarf galaxies. Their starbursts may 
have been triggered by past interactions with their companions (M81 for NGC3077, and 
NGC5236 for NGC5253), but these interactions occurred a few hundred Myr ago, and 
are unlikely to be the direct cause of the starbursts observed today. 



3. Charting the Shock Ionized Gas 

The stellar population study of these galaxies was complemented with the investigation 
of the conditions of the ISM within and around the four starbursts. The standard emission 
line diagnostic diagram log([OIII]/H/3)-log([SII]/Ha) was used to discriminate photoion- 
ized from non-photoionized gas ((Baldwin, Phillips fc Terlevich 1981| ). HST/WFPC2 im- 
ages were obtained with narrow-band filters in correspondence of the forbidden lines 
[OIII](5007 A) and [SII](6726 A), and of the hydrogen recombination lines Ha and H/3. 
The diagnostic diagram was constructed by dividing each image in a grid of bins 3-6 pc 
in size, measuring line fluxes in each bin 5 a above the noise level, and deriving line ratios 
from spatially corresponding bins. Figure [2] (left panel) gives an example of the resulting 
diagrams i C alzetti et al. 2004"j> . 

The 'maximum starburst line' defined by Kewley et al. (2001) was used as the separat- 
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Figure 1. (Left). The color-color diagram, UV— V versus V— I (in HST magnitudes), of the 
stellar clusters in NGC3077. The small open points are the observed data, while the filled points 
are the reddening-corrected data, with the direction of the reddening vector shown as a diagonal 
line at the bottom-left of the plot. The two sets of points are connected by lines. The sizes of the 
circles concentric to the filled points are proportional to the total V-band flux of eac h clusters. 
The colors of an instantaneous burst model from Starburst99 jLeitherer et al. (1999)} are shown 
for the age range 1-900 Myr as the continuous line across the diagram. (K.ight). The age-mass 
plot of the stellar clusters in NGC3077. Each data point is shown with its e rror bar. The curve 
at the bottom represents the 90% completeness limit. Both figures are from Harris et al. 20041 

ing line between photoionized gas (to the bottom-left of the line) from non-photoionized 
gas (to the top-right of the line, Figure^- This can be considered a conservative def- 
inition of non-photoionized gas, yet it yielded clear identifications in all four galaxies. 
The non-photoionized gas in these galaxies shows two basic morphologies: enclosed bub- 
bles and shells/filaments (Figure right panel), and is likely due to shock ionization 
( Martin 1997). Indeed, the presence of star-formation-induced shocks is supported by ob- 
servations of diffuse X-ray emission in the centers of all four galaxies. The shock-ionized 
component is responsible for a small fraction, about 3%-4%, of the total Ha luminosity 
(Tabled column 7), but is spread over a significant fraction of the area covered by the 
Ha emission, about 15%-25%. 

The availability of multi-color (UV-to-I) and Ha images allows us to place contraints 
on the characteristics of the ionizing stellar population, and, once coupled with models, 
on the expected amount of mechanical energy coming out of the starburst. The Ha 
luminosity predicted from the mechanical starburst output can then be compared with 
the observed Ha luminosity associated with shocks. The last column of Table^shows the 
predicted-to-observed Ha luminosities from shocks, with the range marking starbursts 
of different durations, from 10 Myr (left-hand-side figure) to 100 Myr (right-hand side 
figure). In this column, a figure smaller then 1 means that the starburst does not produce 
enough mechanical energy to account for the observed Ha luminosity of shocks, while a 
figure larger than 1 means the starburst can support the observed shock luminosity. 

For the three dwarfs, between 70% and 100% of the mechanical energy is deposited 
within the immediate surroundings of the starburst, roughly within the inner 0.2-1 kpc. 
Furthermore, starbursts with durations shorter than ~30 Myr do not produce enough 
mechanical energy to account for the observed Ha luminosity of shocks. It is worth 
reminding that we adopted a conservative criterion for identifying shocked gas in the 
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Figure 2. (Left). The log([OIII]/H/?)-versus-log([SII]/Ha) diagnostic diagram for NGC3077. 
Each data point represents a bin 0.3" x 0.3" in the WFPC2 images, or about 5.6x5.6 pc 2 at the 
distance of the galaxy (Table 0. Typi cal 1 a uncertainties a re shown at the bottom right of the 
plot. The 'maximum starburst' line of |Kewley et al. (2001)| is shown as a continuous line. This 
line is used to separate areas dominated by photoionizing processes (below-left of the line) from 
non-photoionizing processes. (Right). The bins identified as dominated by non-photoionization 
(shocks) in the diagram at left are reported in black on a gray-scale Ha image of NGC3077. 
North is indicated by the arrow. The position of Ha shells is shown by white poligons, with 
their naming convention (IMartin 19 98 ). while the posit ions of the CO emis sion (Walter et al. 
2002) are shown as white ellipses. Both figures are from lCalzetti et al. 20041 



galaxies, thus the limitation on the minimum duration of the starbursts is quite stringent. 
In addition, the result is independent of whether the parent galaxy is isolated or in 
interaction with other galaxies. Since the vast majority of the starburst mechanical output 
is recovered in the immediate surroundings of the starburst itself, it is not unlikely that 
the shocks can trigger renewed star formation. This inference could help explain the long 
star formation timescales required by the high luminosity of the Hcv emission associated 
with shocks, and the even longer timescales derived from the stellar population studies. 

For the giant spiral, only between 35% and 75% of the available mechanical output is 
recovered in the surrounding ISM, for a variety of starburst conditions. Thus, there are 
no restrictions on the current starburst coming from the characteristics of the shocked 
gas. In addition, the non-photoionized gas is located in two main bubbles, and there 
is a striking absence of extended features (shells or filaments). This suggests a confined 
starburst. 



4. Conclusions 

The characteristics of the stellar populations and ISM in the starbursts of four nearby 
galaxies suggest a strong coupling between the starbursts and the surrounding ISM. 
This is particularly true for the three dwarf galaxies in the sample. Here the amount of 
observed Ha luminosity associated with shocks corresponds to about 70%-100% of the 
total mechanical output from the starbursts, and places a lower limit of about 30 Myr 
to the duration of the starbursts (i.e., these cannot be 'bursts'). Such large fractions of 
mechanical energy available in the immediate surroundings of starbursts may suggest 
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triggered star formation as a mechanism to explain the long duration timcscales inferred 
from stellar population studies («100-300 Myr), in the absence of direct, recent triggers. 
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